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Abstract 


This is the second semi-annual report for NAS5-97039. Work in this project is 
related to NAS 1-20666, also funded by NASA ACMAP. The overall objective of this 
project is to improve the understanding of coupling processes between atmospheric 
chemistry and climate. Model predictions of the future distributions of trace gases in the 
atmosphere constitute an important component of the input necessary for quantitative 
assessments of global change. We will concentrate on the changes in ozone and 
stratospheric sulfate aerosol, with emphasis on how ozone in the lower stratosphere would 
respond to natural or anthropogenic changes. 

The key modeling tools for this work are the AER two-dimensional chemistry- 
transport model, the AER two-dimensional stratospheric sulfate model, and the AER three- 
wave interactive model with full chemistry. We will continue developing our three-wave 
model so that we can help NASA determine the strength and weakness of the next 
generation assessment models. 
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(I) Two-dimensional Chemistry-Transport Model 


(1. 1) Model Enhancements 

The following model enhancements were implemented in the past 6 months, 
primarily for the purpose of performing the WMO/UNEP ozone trends calculations. (1) 
We downloaded from the Goddard DAAC reanalyzed temperature data from NCEP for the 
1987-1995 period. This complements the 1979-1986 data we had obtained in 1996. The 
data was processed into monthly average zonal means and monthly temperature probability 
statistics on the AER 2-D CTM grid. We have used the data both to generate climatological 
temperature and temperature distributions for the entire 1979-1995 period, and for 
calculation of historical ozone trends employing the actual temperature for each year 

between 1979 and 1995. (2) The reaction CIO+OH— >HCI+0 2 was added to the 2-D CTM. 

Though it has been know for some time that use of this reaction would improve agreement 
between models and measurements near 40 km, a definitive measurement of the rate 
constant was not made until 1997 when Lipson et al. (1997) measured the branching ratio 
for the CIO+OD reaction and found 5-6% branching to DC1. (3) Five additional CFC and 
HCFC species were added to the model in order to account for future CFC substitutes in 
our ozone trends calculation. These were CFC-1 14, CFC-1 15, HCFC-141b, HCFC- 
142b, and HCFC-123. 

(1.2) Studies and Publications related to the 2-D CTM 

(1.2.1) ODP ofCH ^r 

We submitted a paper to the Journal of Geophysical Research on the ozone 
depletion potential (ODP) of CH 3 Br. The ODP of CH 3 Br is determined by combining the 
model calculated ODP/BLP ratio (where BLP is the bromine loading potential) for CH 3 Br 
and its atmospheric lifetime. This paper examines how changes in several key kinetic data 
affect the ODP/BLP ratio. The key reactions highlighted in this study include the reaction 
of BrO+H0 2 , the absorption cross-section of HOB r, the absorption cross-section and the 
photolysis products of BrON0 2 , and the heterogeneous conversion of BrON0 2 to HOBr 
and HN0 3 on aerosol particles. By combining the calculated ODP/BLP ratio with the latest 
estimate of 0.7 year for the atmospheric lifetime of CH 3 Br, the likely value of ODP for 
CH 3 Br is 0.39. The model calculated concentration of HBr (~ 0.3 pptv) in the lower 
stratosphere is substantially smaller than the reported measured value of about 1.5 pptv. 

The model can reproduce the measured value of 1.5 pptv if one assumes a yield for HBr of 
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1 .3% from the reaction of BrO+OH, or a yield of 0.6% from the reaction of BrO + H0 2 . 
Our calculations show that the effect of these assumed rates on the model calculated 
ODP/BLP ratio is minimal: practically no impact for the assumed BrO+OH yield and 10% 
smaller for the BrO+H0 2 case. 

A copy of the manuscript is included as an appendix. 

(1.2.2). Adjusting Kyy in 2-D CTM according to observations 
a) Using exchanging rates to adjust in subtropics 

The subtropical barrier reduces the mixing between the tropics and the mid-latitudes 
in the lower stratosphere, and has significant effects on tracer distribution in the 
stratosphere. It is important for 2-D CTMs to have a realistic exchange rate between the 
tropics and mid-latitudes to simulate the effects of the subtropical barrier. The exchange 
rate is determined by the eddy diffusion coefficient K yy in the model. In an earlier study, 

we introduced a tropical pipe by reducing the K yy in the tropics from 0.3 to 0.03 (10^ 

m 2 /sec), (Weisenstein, et al. 1996). Recently, in situ measurements of chemical species 
with a wide range of local lifetimes have been used to quantify the air exchange rate 
between the tropics and mid-latitudes in the lower stratosphere from the tropopause to 
around 21 km. It is found that the mid-latitude air is entrained into the tropical lower 
stratosphere with a replacement time scale of 10-18 months (Minschwaner et al., 1996; 
Volk et al., 1996). Meanwhile, Schoeberl et al. (1997) estimate the mixing time between 
the tropics and mid-latitudes in the 20-28 km region to be at least 18 months using the QBO 
signals in N 2 0/CH 4 ratio and tropical winds from UARS measurements. These results 
demonstrate that the tropical barrier allows a moderate penetration. 

The horizontal diffusion fluxes and the mixing ratio of CCI4, CF2CI2, CFCI3, 
CH3CCI3, CH4, HNO3, N2O, NOy, and O3 calculated in the AER 2-D model are used to 
derive the exchange rates between the tropical and the mid-latitudes in the lower 
stratosphere. They are compared with the exchange rates estimated from observations. 

Using the original model (global diffuser model) with a K yy of 0.3 x 10 6 m 2 /sec in the 
lower tropical stratosphere the model calculated exchange time scales are around 5-6 
months. The tropical pipe model results show a much slower exchange rate of 38-60 

months. Our studies showed that values for K yy of 0.13 x 10 6 m 2 /sec from the tropopause 


Atmospheric and Environmental Research, Inc. 
NAS5-97039 Semi-Annual Report 


3 


to 21 km and 0.07 x 10 6 m 2 /sec from 21 km to 35 km in the model (leaky pipe model) will 
generate an exchange time scale of around 14 months below 20 km and 18 months in the 
20-30 km region, in good agreement with the estimates from observations. 

The model calculated lifetimes for species like N 2 O and CFCs in the leaky pipe 
model are reduced by about 10% compared to the global diffuser model. Our model results 
show that the calculated ozone response to HSCT engine emissions is also smaller in the 
leaky pipe model. 

The paper that summarizes the results of this study was sent to Journal of 
Geophysical Research for publication. Two reviewers have made some suggestions and 
we are in the process of revision. The revised version was sent back to Journal of 
Geophysical Research in February 1998. 


b) Using latitudinal gradient of N 2 0 and NOy/0 3 to adjust Ky y near polar vortex 

It is well know that the isolation of the polar vortex is one of the preconditions for 
the Antarctic ozone hole to be developed (WMO 1988). The mixing ratio of long-lived 
tracers like CH4 and N20 usually are quite different inside and outside of the vortex due to 
the isolation of the vortex and the descending motion of the air inside the vortex during the 
polar winter night. The potential vorticity (PV), which is a good dynamic tracer, has been 
used to define the boundary of the vortex where the PV has large horizontal gradients 
(Manning et al., 1995). The gradient of the long-lived tracers across the vortex is usually 
underestimated in 2-D model simulations. There are two dynamical parameters in AER 2-D 
CTM which affect the horizontal gradient of long-lived species around the vortex, K and 
the circulation. We found that the value of K yy currently used in AER 2-D CTM during the 
winter high latitudes is probably too large to maintain sharp gradients for long-lived species 
near the polar vortex. We put a polar barrier in the model by reducing Kyy from 10 10 
cm 2 /sec to 10 8 cm 2 /sec between 50° to 70° during the winter. It is found that the horizontal 
gradient of N 2 0 and the ratio of NOy/0 3 near the vortex is quite sensitive to the change of 
Kyy. Fig. 1 shows the model calculated annually averaged ratio of NOy/0 3 at 64 mb with 
and without the polar barrier along with the observations. The polar barrier seems to have 
improved the model simulation. The generalized coordinate method, reported at the M&MII 
workshop (1997) by David Considine of GFSC, used to process UARS N 2 0 data 
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produces larger gradient of N 2 0 near the polar vortex than ordinary zonally average does. 
We will use this method to process UARS N 2 0, NO y and 0 3 data and use their gradient 
near the polar vortex to fine-tune K yy in the polar region. 


(II) 2-D Interactive Model 

(II. 1) Improvements of the 3-wave Model 

The current version of the AER 3-wave interactive model has a vertical resolution of 
3.5 km. In order to meet the future requirement for assessment studies, the resolution 
should be increased especially near the tropopause so that the exchange between the 
troposphere and the stratosphere can be better simulated. We decided to increases the 
vertical resolution of the model by three times to match the fine resolution of 1 .2 km used 
in the AER 2-D CTM. The model with finer resolution is running on a SGI machine. The 
results show some differences in the temperature fields between two versions of the model 
with different vertical resolution. The causes of the differences could be the scale 
dependent parameters in the dynamics module such as the sub-grid fourth order diffusion 
coefficients. We are adjusting these parameters to improve the model with finer vertical 
resolution. 

The quasi-biennial oscillation (QBO) is the primary mode of variability in the 
tropical lower stratosphere and has been well documented (e. g. Reed et al., 1961). The 
QBO signals in the temperature have also been observed in the same area (e. g. Reed, 

1964). Associated with the tropical temperature QBO anomalies is a meridianal circulation 
in the tropics and subtropics, which should have a direct effect on the tracer distribution in 
the lower tropic stratosphere. QBO signals in total ozone have been observed for over 
three decades (Funk and Gamham, 1962). The distributions of N 2 0 and CH 4 observed by 
UARS/HALOE also show QBO signals (Schoeberl, et al., 1997). In order to improve our 
model’s ability to analyze the observed data or to use the model to assess atmospheric 
effects of pollutant, QBO should be included in the model. It is difficult for our 2-D 
interactive model to generate tropical waves to provide correct momentum sources for QBO 
due to the coarse vertical resolution and the simplified tropospheric dynamics. So instead, 
the zonal-mean QBO winds, averaged from NCEP data, are imposed on the model by 
adding a relaxation term in the zonal-mean zonal momentum equation. The time scale of 
relaxation is about 100 days. The model generates the right QBO in the zonal winds but 
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only a weak QBO signal in the temperature. Accordingly, the QBO signals in ozone and 
other tracers are small. We are analyzing the model outputs to diagnose the causes of the 
weak QBO signals in temperature. 

(II. 2) Coupling of the 3-wave Model to the 2-D 19-Zone Climate Model 

In order to examine the total effects of solar variation and stratospheric aerosols, 
including their indirect effects from changes in ozone, we are upgrading the 2-D climate- 
chemistry model (CCM). One of the refinements is to replace the 2-D dynamic-chemistry- 
transport module of the 2-D CCM with the 3-wave model (Shia et al, 1997), which 
resolves the three longest planetary waves and calculates the eddy diffusion coefficient for 
chemical tracers from the resolved planetary waves. Thus, the model includes more 
interactions between chemistry, radiation and dynamics. In addition, the 3-wave model has 
updated heterogeneous chemistry. 

The new model is running well on a multi-processor SGI machine when the data 
exchange between the two modules, the 3-wave model and the climate model, is turned off. 
The next step is to establish full interactions between these two modules 

(III) Assessment Related Activities 
(ffl.l) Models and Measurements Workshop (II) 

Our participation in the Models and Measurements Workshop has involved several 
of us at AER in different ways. Malcolm Ko contributed to the definition of one of the 
experiments. Courtney Scott is responsible for the analysis of the source gases 
experiment. The group, as a whole, has produced model results for the UADP archive. 
Malcolm Ko attended the first meeting in August and both Malcolm and Courtney attended 
the second meeting in November, 1997. At the meeting, Courtney made a presentation of 
the results from the source gases experiment. A write-up of the experiment was also 
submitted to UADP in January, 1998. 

Malcolm Ko will also help co-ordinate the preparation of the final report. 

(in.2) Ozone Assessment Report 
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Malcolm Ko is a co-author for Chapter 1 of the WMO Ozone Assessment Report 
with lead authors for the chapter, Ronald Prinn and Rudolf Zander. Chapter 1 deals with 
long-lived trace gases. Malcolm Ko worked with Michael Volk to prepare a section on the 
atmospheric lifetimes of the gases as can be determined from model calculations and 
observations. The model calculated lifetimes are taken from the Models and Measurements 
Workshop (II) calculations. The results from observations are taken from Volk et al. 

(JGR, 102, 25543, 1997) which uses historical records of tropospheric measurement to 
correct the observed correlation in the lower stratosphere to obtain steady state lifetimes. 
One result of the study is that there are indications that the lifetime of CFC-1 1 may be 
shorter than the currently accepted value of 50 years. A clearer picture should emerge 
when the uncertainty in the lifetime can be reduced. 

Jose Rodriguez is lead author for Chapter 2 along with Michael Kurylo. Jose’s 
effort is funded by this project. Chapter 2 deals with source gases that are removed by OH 
in the troposphere and a large part of the chapter was devoted to methyl bromine (CH 3 Br), 
including its ozone depletion potential 

We also performed some chlorine loading calculations and communicated our 
results to Guus Velder who is one of the lead authors of Chapter 1 1. Finally, results from 
a 1970-2050 trend calculation with the AER 2-D CTM were submitted to Chapter 12 by the 
January 1998 deadline. This calculation was performed with observed monthly mean 
temperatures and temperature probability distributions between 1979 and 1995. 

(IV) Chlorine loading and stratospheric cooling 

Recent winter/spring observations in the Arctic [Manney et al., 1996,1997; Muller 
et al., 1997; Rex et al., 1997] reveal further reductions of the ozone content despite nearly 
complete disappearance of Pinatubo aerosol from the stratosphere [Thomason et al., 1997] 
and essentially constant levels of chlorine and bromine loading [WMO, 1995], Existing 
photochemical models of the stratosphere with heterogeneous chemistry predict little or no 
further stratospheric ozone depletion, unless an alternative mechanism is invoked. One 
such mechanism could be the cooling of the lower stratosphere [Oort and Liu, 1993; 

Randel and Cobb, 1994; Ramaswamy et al., 1996] leading to triggering of polar 
stratospheric cloud (PSC) formation. As a result, chlorine activation by enhanced PSCs 
provides further ozone decline during the winter/spring period. Here, we present sensitivity 
study results of the AER box model for typical conditions in the lower stratosphere at 70°N 
during the December 1 - April 1 period for an idealized parcel with different assumed 
stratospheric coolings and chlorine loadings. Our calculations show that the imposed 
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stratospheric cooling could further deplete ozone and retard its expected recovery even with 
the projected chlorine loading decrease. We show that, for the present conditions, a 1 K 
cooling could provide the same local ozone depletion as an increase of chlorine by 0.40- 
0.95 ppbv for the scenarios considered. Thus, sustained stratospheric cooling could further 
reduce Arctic ozone content and delay the anticipated ozone recovery in the Northern 
hemisphere even with the realization of the Montreal Protocol and its Amendments. 


(V) The POLARIS Campaign 

Malcolm Ko was in Fairbanks, Alaska for the deployment of the Polaris Campaign 
in July. Jose Rodriguez joined the September deployment, using the opportunity to meet 
with Michael Kurylo to work on Chapter 2 of the WMO Ozone Assessment Report. 

One focus of the analysis continues to be the NO x /NO y ratio. The observed 
NO x /NO y is larger than that calculated using a constrained photochemical model. R.-S. 
Gao (NOAA) and Ross Salawitch (JPL) have tried to use a rate constant for OH + N0 2 that 
is 50% slower than the JPL recommendation. We used the AER 2-D CTM to perform 
some calculations to see its effect on model calculated ozone trends. 

A copy of the manuscript is included as an appendix. 

(VI) Anticipated Activities in the Next Six Months 

Activities for the next six months will be dictated by the various assessment related 
activities of both the Models and Measurements Workshop, and of the WMO report. We 
hope to find time to revise some of our manuscripts for publications. 
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Abstract 
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Tl.e ozone depletion pmcntial (OOP, of methyl bromide (CHjBr, can be determined by 
combining die mode, calculated ODP/BLP ratio (where BLP is the bromine ,oad.ng potential) 

or CH 3 Br and „s atmospheric lifetime. This paper examines how changes in several hey 

kinetic data affect the ODP/RT P l- 

ODP/BLP ratio. The key reactions highlighted in this study include the 

reaction of BrO + H0 2 , the absorption cross-section of HOBr. the absorption cross-sect, on and 

■h= Photolysis products of BrONO,. and the heterogeneous conversion of BrONO, to HOBr 

and HNOy on aerosol particles. By combining the calculated ODP/BLP ratio with the latest 

estimate of 0.7 year for the atmospheric lifetime of CH 3 Br. the likely value of OOP for CH 3 Br 

' S °' 39 ' ThC m0de ‘ Ca ' Cl " ated C0 "“" trali °" °f HBr (- 0.3 pptv) in the lower stratosphere is 
substantially smaller than the reported measured value of about I pptv. The model 


can 


reproduce the measured value of 1 .5 pptv if one assumes a yield for HBr of 1 ,3% from the 

ow that the effect of these assumed rates on the model calculated ODP/BLP ratio is minimal: 
practically no impact for the assumed BrO+OH yield 


and 10% smaller for the BrO+HO 


2 case. 



n.2 


Introduction 


t measurements of the organic bromine compounds in the upper troposphere and 
lower stratosphere have confirmed .ha, CHjBr is the major contributor ,o stratospheric bromine 
mg (see e.g. Wamsley et al., 1998). The ozone depletion potential (ODP) of CH,Br 
provides a measure of the expected effect on stratospheric ozone from emrssion of CH jB r to the 
atmosphere. The ODP of CH.Br is a function of the number of bromine rad.cals de.ivered to 
the stratosphere per uni, mass of CH.Br emitted a, the ground and the ozone removal efficiency 
of the released bromine radicals. Analysis of the former is complicated by the large 
uncertainties associated with the atmospheric lifetime of CH,Br. Apart from its removal by 
reaction with OH in the atmosphete, CHjBr is also removed by deposition to soil and oceans 
Aithough the impact of bronune chemistty on stratospheric ozone has been studied for close to 
two decades (see review in WMO. 1994), our knowledge of atmospheric bromine is still 
limited by the lack of direct observations of some of the bromine 
With BrO being the only bromine species observed 


species in the stratosphere, 
in aircraft campaigns, the type of process 
Studies that can be used to verify our understanding of the partitioning in the atmosphere is 

limned. The recently reported measured HBr concentration of about 1 pptv (Nol. e, al„ 1997; 

Carlotti e, al„ 1995, Johnson e, al., 1995; Park e, al.. 1989) in the lower stratosphere further ’ 

complicated the picture as models are no, able to reproduce this value with cutren, rate data 
recommendations. 


Changes in the values of the ODP of CH jBr in the pas, have resulted from changes in 

both reaction rate constants for bromine radicals and changes in the atmospheric lifetime of 

CHjBr. In this paper, we review how the calculated ODP for CH,Br has changed with 

mphasts on the recent changes due to the revised rate constants in JPL-97 (DeMorc 
et al., 1997). 
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2. OOP for CH 3 Br 


The Ozone Depletion Potential (ODP), first introduced by Wuebbles (1983), is now 
extensively used by policy makers as a simple relative measure to assess the potential impact of 
a given compound on the stratospheric ozone layer. The ODP is calculated relative to CFC-1 1 
and depends principally on two major factors: a)the bromine loading potential (BLP) which 
measures the amount of bromine radicals delivered to the stratosphere per unit mass emitted at 

the ground, and b) the ODP/BLP ratio which measures the relative catalytic efficiency of the 
radicals in stratospheric ozone removal. 


The bromine loading potential (BLP) reflects the amount of inorganic bromine added to 
the atmosphere from a given emission. An estimate for the amount of inorganic bromine 
present in the atmosphere due to the emission of a particular bromine-containing compound is 
given by the tropospheric concentration of the compound multiplied by the number of bromine 
atoms per molecule. At steady state (i.e., when additions to the atmosphere equal removal 

from the atmosphere), the estimate for the bromine concentration in parts per billion by 
volume (ppbv) is given by 


[B r y]sSx — E SSx T x Nge A x 


where N Br is the number of halogen atoms per molecule of compound X, T x is the atmospheric 
lifetime of compound x, E SSx is the steady state emissions in mass per unit time and A. is: 


Hx mass of the atmosphere lu 

where p is molecular weight in grams, and the mass of the atmosphere is 4.8x 10 21 gm. The 
BLP is the ratio of the steady state concentration of bromine released by compound x for a 
given emission rate to the steady-state concentration of chlorine released by CFC- 1 1 for an 
equal (by mass) emission rate: 


BLP = 



The ODP is defined as the ratio of the ozone depletion produced by a unit-mass 
emission of compound X to the ozone depletion produced by a unit-mass emission of CFC- 1 1 
(the traditional reference gas). Numerical models that incorporate the pertinent chemistry and 
trace gas transport of the atmosphere are used to compute the ODP. For such modeling 
exercises, the steady state ODP for compound X can be given by (Fisher et a]., 1990): 


calculated steady state ozone depletion of about 1% due to X 

emission rate of X to produce this decrease of ozone in the mod el 
calculated steady state ozone de pletion pfaboui 1% due To CFC - 1 1 — 

emission rate of CFC- [ 1 to produce this decrease ofozone in the mode/ 
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where 1% depletion is required to yield results above model noise levels but with,,, the region 


of linear model response. The ratio BEF X 


defined by 


ODP~ 

ls a measure of the relative 
-* X 


ozone removal efficiency of the bromine radicals delivered to the stratosphere. In the 
definition, the behavior of CFC-1 1 is used as the reference standard. The value of BEF 
corresponds to the spatial and seasonal average of the local ozone removal efficiency of the 
bromine radicals released by X relative to the chlorine atoms released by CFC-1 1 . Note that 
the BEF value is species dependent, i. e„ if two species release bromine radicals at different 
locations in the lower stratosphere, the BEF values will be different for the two species. 


Because of differences in model formulation, different models will calculate slightly 
different atmospheric lifetimes. This is further complicated by the fact that some species are 
affected by surface removal processes that are difficult to model. Previous analyses (see e.g. 
Solomon et al„ 1992) have shown that the model calculated ODP/BLP ratio depends on the 
local stratospheric lifetime of the source gas. Changes in lifetimes due to tropospheric and 
surface processes would not affect the ODP/BLP ratio as long as they are not so fast that the 
mixing ratio of the source gas becomes non-uniform in the troposphere. Making use of this 
property, one can define the ODP and BLP using the standard recommended lifetime ( Trec ) as 

follows. The model is used to calculate atmospheric lifetime (x modeI ) and (ODP) model 

according to its definition, taking care that all the stratospheric processes are properly 
accounted for. A value for (BLP) model is then defined as 


(BLP) mo deI = 


m 


£ x modelJx 


I. p ] CFC-i i f T model]cFC- 


and (he model calculated BEF X j s g ; vcn by 


BEF X = ^P Dl ^moAd 

<**■0.0 drf 


The ODP value based on the recommended lifetime is 


defined as 


(ODP)xrcc = (BLP) Trcc 

<**A»od«f 


- (BLP) Trec BEF X 


where (BLP) Xrcc j s defined as 


(BLP) Trec = 




and 


Tree is determined through consensus of the modeling 


community. 


Table , summarizes how lh e calcu.ated va,ues of ODP for CH,Br has changed 
b- These values are calculated by the 


since 


AER model using different rate 


constants as available 


at the time. Between 1990 and 1997, the BEF for CH R h - 

CH 3 Br has increased while its lifetime has 

eercased. On balance, the ODP has decreased from 0.99 to 0.39. 
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Atmospheric Lifetime of Methyl Bromide 

The net atmospheric lifetime (t, oc ) for methyl bromide is given by: 


1 


T 


101 


1 


gas 



T T 

ocean ^land 


the above expression, x gas denotes the lifetime against atmospheric removal by gas-phase 

reactions and photolysis of CH 3 Br, and is defined by the total column density (cm-*, divided by 
the integrated loss rate (cm V*) over the whole atmosphere. The lifetime against ocean 
uptake, Wan, , s defined as the methyl bromide atmospheric burden divided by the flux into the 

ocean. A similar definition holds for the lifetime against uptake by land surfaces, T Iand . 


Although photolysis and reactions of CH 3 Br with O('D), 0(3 P)> N o 3 , and Cl take place 
in the atmosphere, gas-phase removal of methyl bromide is likely to be dominated by its 

tion with tropospheric OH. Based on the latest recommended rate constants of the reaction 


of CH 3 Br with OH, and of CH 3 CC1 3 with OH, together with the most recent analysis of 
ALE/GAGE data (Prinn et aL, 1995), a partial lifetime due to reaction with OH in the 


troposphere for CH 3 Br of 2. 1 ± 0.5 years, is derived. 


When combined with the partial lifetime 


due to stratospheric removal, the estimated T eu is 1.7 years. 


Hydrolysis and reaction with chloride ions in sea water constitute an important sink for 
atmospheric CH 3 Br. This is in contrast to halocarbons and proposed substitutes which, due to 

low solubility (typically H - IQ-2 moles/l/a,m)and aqueous degradation rates (typically k w - 


i 


1 


10-9 s-'), have ocean uptake lifetimes of order of hundreds of years (Wine and Chameides, 
1990). Butler (1994) described a method to calculate the (partial) atmospheric lifetime of 

CH 3 B, due to removal by the ocean (T occan ) in terms of the surface mixed layer depth of the 

ocean, sea-surface temperature, wind speed, thermocline diffusivity, solubility of CHjBr in sea 
water, and degradation rate coefficient in sea water. The Butler (1994) estimate was carried out 
using a two-box model so that weighted averages of the parameters were used in the 

calculations. The calculated best estimate for x occan was 3.7 years, with a possible range from 


1.4 year to 14 years. These values were adopted in the WMO 1994 report. The same 
methodology was used in Yvon and Butler (1996) but with the parameters calculated 
dtvidually for each 2 by 2 grid over the oceans. Thus, they were able to take into account 
the co-variations of the parameters before summing the effects to obtain the global partial 

lifetime. With this method, they obtained a best estimate for r oceJ „ of 2.7 years, with a range 

from 2.4 years to 6.5 years. The latest update (Yvon-Lewis and Butler, 1997) includes 

biological removal and lowers the lifetime 1,8 years, with a range from 1.1 yearn to 3.9 
years. 


Measurements of deposition velocity of CH 3 Br over land are not available. Based on 
the uptake rates measured for five types of soil. Shorter e. al. (1995, computed the partial 
lifetime of CH 3 Br. The value is based on the sum of the annual uptake fluxes (g yr I) by the 
five soil types, each of which is estimated as the product of the uptake rate for the soil type (g 


m-2 d-l). the area of the soil type on the globe (m=), and the number of active days in a year (d 
yr'). The five soil types were assigned two uptake rates, corresponding to deposition 


velocities of 0.12 cm s ' and 0.02 cm S I 


The value obtained due- to soil removal is 3.4 years. 


Combmmg the three partial lifetimes provides a best estimate of 0.7 


years. 
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4 Partitioning of the Bromine Radicals 

We used the AER 2-D model (Weisesntein et al., 1996) with the standard JPL reaction 
rate constants for the gas-phase reactions and parameterization for heterogeneous reactions on 
sulfate aerosol. Following the method used in WMO (1992), the ODP is calculated based on a 
background chlorine and bromine concentration of 3.5 ppbv and 19 pptv, respectively. 


The concentrations of the bromine radicals at 38°N for spring conditions calculated using 
JPL-97 recommended rate constants are shown in figure 1. Note that at that latitude, BrONO, 
and BrO are the most abundant species. The concentrations of HOBr and HBr are less than 1 
pptv. Figures 2a and 2b show the computed [BrO ]/[Br y ] ratio along with observations reported 
by Avalonne et al. (1995). It is evident that the BrO concentrations calculated using JPL-97 are 
larger than those calculated using JPL-94. This accounts for the larger BEF in Table I 

The reported measured concentration of 1.5 pptv for HBr between 25 and 35 km 
(Carlotti et al., 1996) is much larger than the model computed value of about 0.3 pptv. In 

section 6, we will discuss ways of increasing the calculated concentration of HBr in the model 
and how they affect BEF. 


5. Sensitivity of BEF to Rate Data 


The following reactions represent the primary bromine-driven catalytic ozone removal 
cycles in the stratosphere: 


CIO + BrO + hv — » Br + Cl + 0% 
Br + O^-* BrO + 0 2 
Cl + 0 3 -4 CIO + 0 2 


(I) 


p 10 


CIO + BrO + hv -> Br + Cl + 0 2 
B t -f- 0 3 — > BrO + O 2 
Cl + 0 3 -> CIO + 0 2 

2O3 -» 30 2 . 


and 


BrO + H0 2 -> HOBr + 0 2 ' 
HOBr + hv -> OH + Br 
Br + O3 — > BrO + 0 2 
OH + 0 3 — » H0 2 + 0 2 , 


20 3 30 2 . 

BrO + O — > Br + 0 2 
Br + O3 — ) BrO + 0 2 


O +O3 — > 20 2 . 

2(Br + 0 3 — > BrO + 0 2 ) 1 
BrO + BrO — > 2Br + 0 2 J 


20 3 30 2 . 


BrN0 3 + hv -> Br + N0 3 
NO3 + hv — > NO + 0 2 
Br + O3 BrO + 0 2 
NO + O3 — > N0 2 + 0 2 
BrO + N0 2 + M ^ BrN0 3 + M 


> 


J 


20 3 -> 30 2 . 


(II) 


(HI) 


(IV) 


(V) 


Current rate constants indicate that, in the present day atmosphere, the first three cycles 
contribute about an equal amount to ozone depletion in the lower stratosphere, with the last two 
about an order of magnitude smaller. 
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To illustrate how changes in rate constants affect BBF, we summarize the key revisions 
between the JPL-94 and JPL-97 in Table 2. Changes to reactions (R I ) through (R8) (see Table 
2) impact the ozone removal efficiency either directly because they are the rate limiting step in 
the cycle, or indirectly by changing the partitioning between BrO and the reservoir species. We 

explore the dependence by performing a series of calculations to obtain the BEF in each case. 
The results are summarized in Table 3. 


Reaction Rl in Table 2 has received much attention. A key revision in JPL-92 that 
affects the BEF value for CH 3 Br is the faster rate constant for reaction of BrO and H0 2 . With 
an estimated temperature dependence, the rate recommended in JPL-92 based on room 
temperature value by Poulet et al. ( 1992) is about a factor of 10 faster at stratospheric 
temperatures than that recommended by JPL (1990). Reaction Rl represents the rate limiting 
steps in cycle (11). With the rate constant from JPL-92 (a ten-fold increase compared to JPL- 
90), the BEF increases from 37 to 49. 

In the JPL-94 calculation, we included the absorption cross-section for HOBr as 
reported by Orlando and Burkholder (1995), which is the first measured absorption cross- 
sectron. Prior to that, the much larger cross-section of HOCI was used as a proxy, resulting in 
computed HOBr concentration of less than 0. 1 pptv at 1 6 km and 0.8 pptv at 30 km. With the 
Orlando and Burkholder cross-section, the concentration of HOBr is around 1 pptv at 16km 
and 3 pptv around 30 km (see the curve labeled JPL-94 in Figure 3), comparable to those of 
BrO and BrONO,. Inclusion of this cross-section lowers the BEF from 49 to 44. 


Several revisions in JPL-97 identified in Table 2 also have impacts on the model 
calculated HOBr. A smaller rate constant for R 1 results in smallec HOBr (see Case B in Figure 
3). The heterogeneous reaction R3 convert BrONO, to HOBr and BrO. As a result of this 
reaction, both HOBr and BrO are larger (see Case D in Figure 3). The combined effect of R I 
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and R3 is represented in Figure 3 as Case E. For comparison, we performed Case F where we 
have updated ail the JPL-97 rate data except the absorption cross-sections. The value obtained 
in the lower stratosphere is similar to Case E. The difference between Case E and Case F 
around 30 km is the result of the slower rate for R8 (O + HOBr) in the JPL-97 
recommendation. Figure 4 shows the absorption cross-section of HOBr from Orlando and 
Burkholder (1995) and the JPL-97 recommendation. As noted in JPL-97, the recommended 
cross-section is influenced by the measurements reported by Rattigan et al. (1996). With the 
larger cross-section between 300 nm and 400 nm and the extension of the cross-section beyond 
400 nm, the photolysis rate of HOBr is three times larger in the lower stratosphere leading to a 
smaller concentration of HOBr (see cases F and JPL-97 in Figure 3). 

We will now discuss the effects of the rate revisions on the model calculated BEF for 
CH,Br. The Orlando and Burkholder cross-section decreases the sensitivity of the model 
calculated BEF to changes in the reaction rate constant for reaction (I). With HOBr being a 
major reservoir, any increase (decrease) in k, is compensated by decrease (increase) in [BrOj 
so that the change in the rate limiting step k , [BrO][HOJ will be smaller than the corresponding 
change m k,. This is illustrated by the results for case B in Table 3. Despite a factor of two 
decrease in the rate constant for Rl. the ODP for case B is practically the same as that for case 
A. Revision of the BrONO, cross-section involves both the magnitude and the product yield. 
The measurement repotted by Burkholder el al. (1995) extends the measured values to beyond 
390 nm. With the new cross-section, the photolysis rate is a factor of 1.25 larger in the lower 

stratosphere. For BrONO, photolysis (R2). only ^ ^ leads to ozone destruction as the 
^BiOtNoi b ™ch is a null cycle. Previously, it has been assumed that <t> =1 w :,k a, 

Br ♦ N03 BrO + N02 

=0. The new values from Nickolaisen and Sander (private communication) indicate that d» 

BrfN03 

0.29 with ^ B(O+N02 = 0.71. Thus the new branching should decrease BEF. However, 
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diagnostics reveals that the cycle (V) contribution is small compared to the other cycles around 
20 km. Thus, the effect of the change in branching on BEF is small as illustrated by Case C. 
However, the expected effect from the extension of the absorption cross-section to longer 
wavelength is to suppress BrON0 2 and increase BrO and BEF (see Case F and Case H). 

Reaction (R3) is a major source of HO, in the lower stratosphere via production of 
HOBr and its subsequent photolysis. The reaction also leads to a higher concentration of BrO 
and smaller concentration of NO,. The smaller NO, provides the coupling to the CIO, and HO, 
cycles. Cases D shows a 19% increase in the BEF. Our diagnostics reveal increased 

contributions from cycles (I) and (II) where a larger reservoir of HOBr and BrO are produced 
via reaction (R3). 

The combined effects of RI, R2b and R3 is given in case E. The resulted BEF can be 
compared with that for Case F, which consists of all JPL-97 updates except the cross-sections. 
The similarrity in BEF for the two cases shows that the effects of R5 through R8 on BEF are 
small, The final four cases in Table 3 illustrate the contributions of the revisions the cross- 
sections on BEF. The HOBr revision has the largest effects among the revisions. The larger 
cross-section leads to larger calculated concentrations of BrOand larger BEF. For BrON0 2> the 
revision to the absorption cross-section again leads to a repartioning of the bromine radicals 
leading to higher concentrations of BrO. The effect from the change in branching is small. The 
cross-section revisions for C10N0 2 and BrCl are each responsible for increasing the BEF by 
0.5 to bring the BEF to 58 in Case J. 
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The ODP and (he BET of CH,Br with HBr concentrations at 1.5 |,,,(v 


The calculated concentration of HBr in Figure I is much smaller than values reported in 
the literature. To obtain calculated concentrations of HBr similar to recent measured values of 

HBr in the lower stratosphere by Carlotli et. al. (1995) and Nolt el al. (1997). we have to invoke 
alternative mechanisms to produce HBr. 


Several studies (Larichev et a!., 1995; Mellouki et al.. 1994) indicated that the product 
s HOBr and O,, although it is not possible to confirm or refute a small yield of HBr by 
direct measurements. Another possibility is a small yield of HBr from the react, on of BrO with 
OH. Two recent publications presented model results to show how a small yield from BrO + 
OH (Chipperfield et al.. 1997, or yields from either reactions (Chartraud and McConnel. 1998) 
could successfully simulate the measured concentration. In two separate cases, we assume a 
1 .3% yield of HBr from BrO + OH. or an HBr yield of 0.6% from BrO + HO,. Annual 
averages of HBr between 25 and 35 km are calculated to be 1.3 pptv and 1.5 pptv for the two 
cases respectively (see figure 5). In both cases, HBr is elevated by over 100% from Case A. 

The concentrations of the remaining bromine species are decreased by up to 10%. The effect 
on BEF is small as evident from the small changes in BrO in figure 6. The calculated BEF 

values are 56 (for the BtO + OH case , and 52 (for the BrO + HO,) case, compared to a value of 
58 calculated with no assumed HBr yield. 
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6. Concluding Remarks 

The calculated ODP for CH 3 Br depends not only on stratospheric bromine chemistry, 
but also on the calculated atmospheric lifetime. Some of the processes which define the surface 
removal are not yet well understood. In contrast to the vast amount of data on C10 x species, 
observations of BrO x species are few. Although observations of BrO in the lower stratosphere 
and the tentative upper limits on HBr above 30 km place some constraints on bromine 
chemistry, they do not rule out models which assume a small (5-10%) HBr branching. 
Laboratory studies of reactions (1) and (2) together with selective (but more frequent) 
stratospheric measurement (e.g., simultaneous measurement of BrO and HBr over an extended 

(20-30 km) altitude region) will help accelerate our understanding of stratospheric bromine 
chemistry. 
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Figure I; Model calculir^ri a 

38-N latitudes. The ca.cu, J ^ ^ ^ ** « 

calculations were perfonned using JPL . 97 ra , c 

RgUre 2 " C ° m P aris °" of model calculated values of the BrO/B ■ 

values derived from the AASE n T ratio compared with the 

model values are for 20 ^77^ " ” * ^ «* <'9,5). The 

0 alt,tudl: and are taken from the 2-D CTM ' 

calculated values of BrY and ozone Th i computed using model 

o«ne. The results are soiled by latitudes. 


Figure 2b: Same as figure 2a 


except that the results are sorted by solar 


Figure 3: Model calculated day-time avers 

Table 3 formed ■ ■ e average concentration of HOBr for di 

lor the description of the cases. 


zenith angles. 

fferent cases. See 


Figure 4: Absorption cross-sections of HOBr f rv, 

JPL-97 recommendation. ' ““ Burkhold =r (1995) and the 
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Figure 6: Same as figure 5 


except for BrO. 



Tabic 1: Values of ODP for CII3Br 



BEF 

key revision 

36 

- 

49 

BrO+H0 2 10 times 
faster 

48 

JPL-94 

44 

HOBr cross-section 
from Orlando and 
Burkholder (1995) 

58 

see table 2 



comment 


removal by OH only 


include removal by 
ocean 


include removal by 
ocean 


include removal by 
soil, and biological 
removal in ocean 
water 



Table 2: Summary of rate revisions for key bromine chemistry reactions in JPL-97 


rate 


Rl: BrO+HO, 


R2a: absorption cross- 
section of B 1 ONO 2 , 


JPL-97 

3.4x1 0- l2 exp(5407T) 


J PL-94 


6.2xl0- I2 exp(500/T) 


temperature dependent 
cross-section extended 
to longer wavelengths 



R2b: quantum yield 


R3: heterogeneous temperature dependent 

conversion of BrONO, rate with y as , arge as 
on suiratc 0 8 
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Tabic 3: Model calculated BEF using different input 
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Figure 1 : Model calculated day-time averaged mixing ratio profiles for spring conditions at 
38°N latitudes. The calculations were performed using JPL-97 rate recommendations 
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Figure 2a: Comparison of model calculated values of the BrO/Bry rati 


values derived from the AASE-II 
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Figure 2b: Same as figure 2a except that the results are sorted by solar zenith angles. 
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Figure 3. Mode! calculated day-time average concentration of HOBr for different cases. See 
Table 3 for the description of the cases. 
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Figure 4: Absorption cross-sections of HOBr from Orlando and Burkholder (1995) and the 
JPL-97 recommendation. 
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Abstract 

Recent winter/spring observations in the Arctic [Manney et ah, 1996,1997; Muller et ah, 

1997; Rex et ah, 1997] reveal further reductions of the ozone content despite near complete dis- 
appearance of Pinatubo aerosol from the stratosphere [Thomason et ah, 1997] and essentially 
constant levels of chlorine and bromine loading [WMO, 1995]. Existing photochemical mod- 
els of the stratosphere with heterogeneous chemistry predict little or no further stratospheric 
ozone depletion, unless an alternative mechanism is invoked. One such mechanism could be the 
cooling of the lower stratosphere [Oort and Liu, 1993; Randel and Cobb, 1994; Ramaswamy 
et ah, 1996] leading to triggering of polar stratospheric cloud (PSC) formation. As a result, 
chlorine activation by enhanced PSCs provides further ozone decline during the winter/spring 
period. Here, we present sensitivity study results of the AER box model for typical conditions 
in the lower stratosphere at 70°N during the December 1 - April 1 period for an idealized parcel 
with different assumed stratospheric coolings and chlorine loadings. Our calculations show that 
the imposed stratospheric cooling could further deplete ozone and retard its expected recovery 
even with the projected chlorine loading decrease. We show that, for the present conditions, a 
1 K cooling could provide the same local ozone depletion as an increase of chlorine by 0.40-0.95 
PPbv for the scenarios considered. Thus, sustained stratospheric cooling could further reduce 
Arctic ozone content and delay the anticipated ozone recovery in the Northern hemisphere even 
with the realization of the Montreal Protocol and its Amendments. 

Introduction 

It is now widely accepted that the spring ozone hole, first detected in Antarctica [Farman 
et ah, 1985; Chubachi et ah, 1985], can be explained primarily due to activation of chlorine- 
containing molecules of anthropogenic origin [Molina and Rowland, 1974] via heterogeneous 
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reactions on PSC surfaces and their subsequent photolysis [Solomon et al., 1986; McEIroy et 
ah, 1986; Molina and Molina, 1987]. The Montreal Protocol with its Amendments are aimed 
to curb the emissions of chlorofluorocarbons (CFCs), eliminating their production after the 
year 2000 [WMO, 1995] thus avoiding further depletion of the global ozone layer [Prather et 
al., 1996]. However, concerns still remain whether an ozone depiction threat can be completely 
avoided. Recent satellite [Manncy et al., 1996,1997; Newman et al., 1997], ozonesondcs, and 
ground-based [Bojkov et al., 1995; Fioletov et al., 1997] measurements clearly show an alarm- 
ing level of stratospheric ozone depletion (~30-40% below its climatological values at 20 kin 
altitude and less than 250 D.U. in total column (i.e. ~30-40% below climatological values)) 
in the Northern Hemisphere during last two winters/springs. According to the current under- 
standing, these large values of ozone depletion are surprising, since the Pinatubo aerosol has 
been mainly removed [Thomason et al., 1997] and the level of chlorine and bromine loading 
remains almost constant [WMO, 1995]. However, more detailed meteorological analysis reveals 
low temperatures in the stratosphere [Manney et al., 1997; Fioletov et al., 1997] and a stable 
longer-lasting polar vortex [Coy et al., 1997] during the past winter. 

The purpose of this study is to investigate the ozone behavior within an idealized typical 
air parcel in the Arctic stratosphere under projected declining chlorine loadings superimposed 
with various assumed cooling trends. We will also show that the time constant associated with 
ozone recovery in the colder stratosphere could be significantly prolonged. We are address- 
ing neither the reasons or dynamical aspects of stratospheric cooling nor the contributions of 
different radiatively-active gases to the stratospheric temperature regime. Also, we are not 
attempting to predict the ozone content above particular locations for particular dates. Our 
sensitivity study tries to capture a typical year-to-year pattern of the 0 3 behavior in the Arctic 
winter/spring stratosphere. We focus on evaluating the chemical implications of stratospheric 
cooling assuming, for the sake of simplicity, the same air parcel trajectory and ignoring ozone 
radiative [Ramaswamy et al., 1996] and dynamical [Rosenfield et al., 1994] feedbacks. Shindell 
et al [1998], using the GISS 3-D climate model with sophisticated radiative and dynamical 
codes but with simplified chemistry, showed that the chemical mechanism plays the dominant 
role in the polar ozone recovery for the projected chlorine and greenhouse gas loadings, thus 
making our analysis particularly relevant. 


2 



Initialization of the AER Box Model 

Photochemical box models have been applied successfully to analyze different aspects of the 
Antarctic ozone hole [e.g. Schoeberl et al., 1996; Danilin et ah, 1996; Kawa et ah, 1997], For 
this study, we use the AER photochemical box model with the heterogeneous chemistry package 
described in Danilin et ah [1996]. We investigate chemical effects in the lower stratosphere, 
since the ozone profile has its maximum there and its perturbation due to PSC processing is 
also strongest there, thus determining to a large extent the magnitude of total ozone column 
depletion. So our box model calculations in the lower stratosphere could serve as a proxy of 
the ozone column change. We simulate the behavior of an isolated air parcel at a latitude 
of 70°N subsiding with a rate of. 1 km/month from ~24 km on December 1 to ~20 km on 
April 1 [Rosenfield et ah, 1994], The latitude of 70°N for our calculations is chosen because 
this latitude is typically located within the polar vortex according to analysis of zonal fields of 
potential vorticity [Manney et ah, 1994]. We constrain initial concentrations in the air parcel 
using the UARS data in a fashion similar to Carslaw et al [1997a]. Initial mixing ratios of 
O3 (3.6 ppmv), HC1 (1.7 ppmv), H 2 0 (5 ppmv), and CH 4 (1 ppmv) are chosen based on the 
HALOE measurements on December 1 at our starting point (70°N, 32 hPa). Initial values of 
CIONO2 (1.2 ppbv) and HNO3 (11.5 ppbv) are taken from the CLAES measurements [Roche 
et ah, 1994]. We assume 3 ppbv of Cl y (=C1 + CIO + HC1 + C10N0 2 + HOC1 + 2(C1 2 0 2 + 
Cl 2 )), 14 ppbv of NOj, (= HN0 3 + NO + N0 2 + N0 3 + HN0 4 + HONO + 2N 2 0 5 + C10N0 2 ) 
and 20 pptv Br y (= Br + BrO + HOBr + BrON0 2 + BrCl + HBr) consistent with the CLAES 
measurements N 2 0 using the N 2 0-NO y [Loewenstein et ah, 1993], N 2 0-C1,, [Woodbridge et ah, 
1995], and N 2 0-Br y [R.Salawitch, personal communication] correlations. Initial partitioning 
among the rest of the chlorine, bromine, and nitrogen species is taken from the AER 2-D model 
[Weisenstein et ah, 1996]. We also assume decrease of NO,, in the air parcel from 14 ppbv on 
December 1 to 10.2 ppbv in late March to be consistent with MLS HNO3 measurements [Santee 
et ah, 1997]. Initial aerosol surface area density (SAD) is equal to 1.8 /rm 2 /cm 3 with assumed 
log-normal distribution with mean diameter of 0.155 pm and width cr=1.49, consistent with the 
aircraft measurements for the background conditions used in Carslaw et al [1994]. The aerosol 
SAD is kept constant during the NAT runs and is variable during the STS runs due to uptake 
of HN0 3 and H 2 0. The shape of the temperature history chosen in this study (see Figure 1) 
mimics the evolution of the vortex-averaged minimum temperature in the lower stratosphere 
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[Zurek et al., 1996] but is increased uniformly by 5 I< to provide ~30% ozone depletion at the 
point of interest during a model run from December 1 to April 1 . This temperature shift serves 
as a ’tuning’ parameter. We also do not superimpose 6-10 day waves on our temperature history 
(like in [Carslaw et al., 1997a]) to avoid additional complications of the day-to-day variability 
of PSCs, which could mask the main point of our study. 

Type I PSC schemes 

One of the key uncertainties for our analysis is what PSC scheme should be adopted for the 
model runs. The issue of how type I PSCs are formed is still a subject for discussion [Carslaw 
et al., 1994, 1997b; Tabazadeh et al., 1994; Tolbert, 1994], However, the model analysis by 
Sessler et al. [1996] shows that, evaluating the chlorine activation and consequent ozone de- 
pletion, a whole spectrum of type I PSC formation possibilities is embraced by the nitric acid 
trihydrade (HN0 3 «3H 2 0, NAT) and supercooled ternary solution (H 2 SO 4 /HNO 3 /H 2 O, STS) 
schemes. Our box model allows type I PSCs formation via NAT without supersaturation (with 
fixed diameter of 1 pm) or via STS [Tabazadeh et al., 1994], The comparative efficacy of chlo- 
rine activation in these schemes also depends on the assumed level of background sulfate aerosol 
and size of NAT particles. Figure 2 addresses differences between these two schemes from the 
point of view of SAD available for heterogeneous reactions in the NAT and STS schemes as 
well as the reaction probabilities and e-folding times of the principal chlorine activation reaction 
CIONO 2 + HC1 — > Cl 2 + HNO 3 . A very steep increase of the reaction probability and surface 
area with decreasing temperature below 196 K drastically reduces the e-folding time of chlorine 
activation from ~200 hours at 196 K to 1.5-3 hours at 191 K. For the particular parameters of 
this study, the STS scheme is more efficient below 191 K, while a steep increase of the the NAT 
SAD between 195 and 192 K makes the NAT scheme slightly more efficient than STS scheme in 
this temperature range. Figure 2 also shows that it is barely possible to use CIO measurements 
to distinguish the NAT and STS schemes for the air parcels exposed to low temperatures for 
a long time (days or so), since all chlorine is expected to be activated in both schemes. Addi- 
tionally, Figure 2 can demonstrate why the ozone hole has been first observed in the Southern 
(but not in the Northern hemisphere) in the 1980s. Colder temperature in the Antarctic lower 
stratosphere (by 7-8 I< comparing with the Arctic stratosphere) and more stable polar vortex 
greatly facilitates chlorine activation, thus producing stronger ozone depletion for the same or 
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slightly lower level of Cly there compared with that in the Arctic. 


Equivalence of stratospheric cooling and increase of chlorine loading 

To investigate the sensitivity of the ozone depletion to a wide range of chlorine loading and 
stratospheric cooling, we performed model runs for 2 to 3.5 ppbv Cl y and 0 to 5 K stratospheric 
cooling. Effects of stratospheric cooling are investigated by lowering the whole temperature 
pattern by 1 to 5 K while keeping its shape unchanged as shown in Figure 1. The change of 
ozone mixing ratio in the air parcel on April 1 as a percentage of its initial value of 3.6 ppmv 
is plotted as a function of chlorine loading and further stratospheric cooling in Figure 3. This 
figure summarizes our findings and underlines the widely accepted point that ozone content 
decreases with increasing level of Cly (draw any horizontal line in Figure 3). However, this 
figure also clearly shows that possible further stratospheric cooling (due to CO 2 or H 2 0 [IPCC, 
1995; WMO, 1995]) has the effect equivalent to an increase of chlorine loading from the ozone 
balance point of view. 

In order to further elucidate results in Figure 3, we introduce "chlorine-cooling equivalent” 
(CCE) as a measure of the change in chlorine loading equivalent to a 1 K cooling for a given level of 
the ozone depletion in the air parcel under specified ambient conditions. Our results for the NAT 
scenario show, for example, that 30% decrease of the ozone mixing ratio in the air parcel could 
result from either 3.1 ppbv Cly and T 6ase , or 2.15 ppbv C\ y and T^g-IK. In other words, a 1 I< 
cooling produces the same effect in the local ozone content as an increase of Cl y by 0.95 ppbv 
(i.e. the cooling-chlorine equivalent is equal to ~ 0.95 ppbv(Cl y )/K). For the STS scenario, the 
ozone sensitivity to cooling is smaller (CCE«0.4 ppbv(Cl y )/K for the 30% depletion isoline). 
The higher values of the CCE for the NAT runs can be explained by the very steep increase of the 
NAT SAD for the chosen present day case (see Figures 1 and 2). However, further cooling leads 
to the saturation of the NAT SAD due to complete conversion of HN0 3 from gas to solid phase, 
thus lowering the CCE values to several tenth of ppbv(CI y )/K. The similar pattern is observed 
for the STS scenario, thus showing that the CCE concept developed here remains valid for the 
different PSC schemes. The absolute values of the CCE is a non-trivial function of the chosen 
ozone depletion, ambient atmospheric conditions, and chosen PSC scheme. The CCE is a local 
parameter and should not be extrapolated to the global scale. However, the concept of CCE, 
applied in more sophisticated global models, could be useful in regulatory discussions. One can 
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also introduce a carbon dioxide - chlorine loading equivalent or ozone depletion potential (ODP) of 
C0 2 which could help to formulate a future regulatory ozone protection policy translating C0 2 
emissions to CFCs emissions after proper scaling using global model calculations. The WMO 
[1995] Report does not consider C0 2 as an ozone depleting species. 

Our calculations for Antarctic conditions (similar to [Danilin et al., 1996] and not shown 
here) show CCE values ~0.15 ppbv(Cl y )/K there. This implies that possible future cooling 
of the Antarctic stratosphere would have a smaller effect on ozone depletion (which is close 
already to its saturation level [Hofmann et al., 1997]) than similar cooling in the Arctic. 

Discussion and outlook 

As a continuation of our sensitivity analysis, we present calculations of ozone behavior 
along the idealized air parcel considered in the previous sections using anticipated future level 
of chlorine loading and assumed stratospheric cooling. We follow the WMO scenarios for the 
01^ (shown by solid black line in Figure 4) and CH 4 concentrations assuming a fixed level 
of Br y (—20 pptv). Figure 4 shows that, ignoring stratospheric cooling, one indeed should 
anticipate a steady recovery of ozone (green lines) in the lower stratosphere immediately after 
reduction of the Cl y loading peaking at the level of 3.08 ppbv at 24 km (or ~3.7 ppbv at the 
top of the stratosphere) in the year 2000. The red dashed line shows “the road not taken” 
[Prather et al., 1996], assuming the level of chlorine is fixed at 3 ppbv and stratospheric cooling 
with a rate of -0.15 K/year. However, even for the relatively modest constant temperature 
trend of -0.05 K/year [Randel et al., 1994; Ramaswamy et al., 1995], ozone remains below its 
1995-level for the next fifty years in the NAT scenario. For the STS scenario, ozone returns to 
the 1995-content near the year 2030. For the larger temperature trend of -0.1 K/year, further 
delay of the ozone recovery is expected with more than 40% ozone depletion, keeping the ozone 
content in the parcel below its 1995 level during next 50 years in the both scenarios. For the 
more severe cooling rate of -0.15 K/year, a larger ozone depletion is calculated (more than - 
40%) with a longer-lasting impact. For the considered scenarios, the maximum ozone depletion 
could be anticipated during the 2010-2020 period for all temperature trends and both PSC 
scenarios. Thus, our calculations show that the ozone depletion in the Arctic stratosphere may 
be deepened and its recovery delayed by several decades when stratospheric cooling is taken into 
account. Both our results and the 3-D climate model calculations of Shindell et al [1998] show 
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that the currently accepted anticipation of ozone healing in the Arctic following the elimination 
of CFCs emissions may be too optimistic. 

Our results show that the Arctic ozone is very sensitive to the ambient temperature in the 
lower stratosphere. We present estimates of only the chemical effects of stratospheric cooling. 
Stolarski [1997] discussed three possible scenarios of future development of the Arctic ozone, 
depending on what the main reason of the observed stratospheric cooling is: a) greenhouse gases, 
b) ozone, or c) just natural variability. Since the natural variability of the Arctic winter strato- 
sphere is large and relevant data sets are short, it is difficult to determine the cause of the 
observed cooling. Each cause could have important implications, a) Carbon dioxide provides 
the strongest cooling in the lower stratosphere among greenhouse gases (e.g. IPCC [1995], 
Clough and Iacono [1995]). If the cooling and the polar ozone hole propagate upwards in the 
near future, it could be considered a manifestation of greenhouse gas effects due to the vertical 
profiles of the C0 2 cooling rate [e.g. Clough and Iacono, 1995]. In this case, the extent of 
polar stratospheric ozone depletion could become considerably larger and its recovery delayed 
by several decades even with the implementation of the Montreal Protocol. If future polar 
measurements and more sophisticated 3-D calculations support our findings, C0 2 emissions 
should be contained for better ozone protection. 

b) Ozone itself is a very important radiative species and there is a strong coupling between 
radiative, chemical, and dynamical processes in the Arctic stratosphere. In this case, more so- 
phisticated 3-D model calculations with all these feedbacks and a state-of-the-art stratospheric 
heterogeneous chemistry module [e.g. Austin et al., 1992] are required to predict future cooling 
and evolution of the possible Arctic ozone hole. 

c) It is also possible that current stratospheric conditions represent a complicated mixture of 
natural variability and human impact. If this is true, further model prediction of the Arctic 
stratospheric ozone recovery is an extremely complex issue. 

Since the current atmospheric conditions in the Arctic are very close to conditions allowing 
further steep ozone loss via an additional PSC formation, a very careful analysis of the con- 
sequences of current or future levels of emissions, which directly (e.g. aircraft H 2 0 and NO x ) 
or indirectly (C0 2 ) facilitate the formation of PSCs or provide additional sites (aircraft sulfate 
particles) for heterogeneous chlorine activation, is required [Peter et al., 1991; Weisenstein et 
al., 1996; DelNegro et al., 1997], Our previous study shows a possibility that severely depleted 
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ozone in the polar region could be diluted to the mid-latitudes after the final stratospheric 
warming [Sze et ah, 1989]. A more detailed study of the reasons for the lower stratospheric 
cooling with particular emphasis on anthropogenic contributions to this effect and consequent 
effects on ozone at polar and middle latitudes will require a use of a 3-D model with all neces- 
sary feedbacks. 
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Figure captions 

Figure 1 . The temperature history adopted for the box model runs at 70°N: solid black 
line show the present day case; black dashed lines depict the present day case lowered by 1, 2, 3, 
4, and 5 K, the almost horizontal dashed lines correspond to the NAT formation threshold and 

the ice frost point, respectively (their slope is due to a subsidence with the fixed H 2 0 mixing 
ratio of 5 ppmv). 

Figure 2. Top panel: Calculated surface area density of STS (solid red line) and NAT 
(solid blue line) particles in ^m 2 /cm 3 (left vertical axis) and reaction probability (right vertical 
axis) of the C10N0 2 + HC1 > Cl 2 + HN0 3 reaction on STS [Hanson and Ravishankara, 1994] 
(dashed red line) and NAT [Hanson and Ravishankara, 1993] (dashed blue line) as a function 
of temperature. Surface area of NAT was calculated assuming its unimodal distribution with 
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diameter of 1 fim and density 1.62 g/cm 3 . Below 205 K, the model of Tabazadeh et al [1994] 
is used to calculate swelling of aerosol particles due to additional uptake of H 2 0 and HN0 3 . 
Bottom panel: Timescale of the C10N0 2 + HC1 — -> Cl 2 + HN0 3 reaction (defined as r 
= 4/ (7 x SurfaceArea x Speedcro N 02 , in hours) on the STS (red) and NAT (blue) shows an 
effectiveness of chlorine activation as a function of temperature. Black symbols show minimum 
temperature during the model runs from 195.25 K (present day case) to 190.25 K (present day 
case-5K). 

Figure 3. Isolines of ozone change in the air parcel on April 1 (in % from its initial value 
of 3.6 ppmv) as a function of chlorine loading and stratospheric cooling at 20 km and 70°N for 
the NAT (top panel) and STS (bottom panel) scenarios. These values correspond to constant 
Br y =20 pptv and the temperature evolution shown in Figure 1 (i.e. zero cooling correspond to 
the T base line and 5K cooling - to the T fease -5K line shown in Fig.l). 

Figure 4. Ozone depletion as function of year for the projected chlorine and methane 
loading and different rate of stratospheric cooling for the NAT (top panel) and STS (bottom 
panel) scenarios. Bromine loading is assumed to be constant at 20 pptv. The black line in 
the bottom panel shows the adopted decrease of chlorine loading at the starting point of the 
air parcel at 32 hPa. We considered scenarios when the present temperature pattern (Fig.l) 
remains constant (green curves) and uniformly lowered at a rate of 0.05 K/year (blue curves), 
0.1 K/year (cyan curves), and 0.15 K/year (red curves). Ozone change in the parcel for the 
constant level of chlorine loading at 3 ppbv for the adopted cooling rate of 0.15 K/year is shown 
by the pink dashed lines. 
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